Raman scattering in a two-dimensional electron gas: Boltzmann equation approach 
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The inelastic light scattering in a 2-d electron gas is studied theoretically using the Boltzmann 
equation techniques. Electron-hole excitations produce the Raman spectrum essentially different 
from the one predicted for the 3-d case. In the clean limit it has the form of a strong non-symmetric 
resonance due to the square root singularity at the electron-hole frequency lu = vk while in the 
opposite dirty limit the usual Lorentzian shape of the cross section is reestablished. The effects of 
electromagnetic field are considered self-consistently and the contribution from collective plasmon 
modes is found. It is shown that unlike 3-d metals where plasmon excitations are unobservable 
(because of very large required transfered frequencies), the two-dimensional electron system gives 
rise to a low-frequency (uj oc k 1 ^ 2 ) plasmon peak. A measurement of the width of this peak can 
provide data on the magnitude of the electron scattering rate. 
PACS: 73.50.-h, 78.30.-j 



Raman scattering is a powerful method for experimen- 
tal studies of elementary excitations in various structures. 
In particular, high-T c superconductors produce Raman 
spectra which remain mysterious over a broad reg-ij 
of frequency. Namely, the high-frequency continuurr 
2A-peak (see, e.g., Ref.El), two-magnon spectra, reveal- 
ing strong mutual influence between antiferromagnetism 
and superconductivityu, still does not have a robust self- 
consistent theoretical description. Therefore, the devel- 
opment of the theory of Raman scattering from different 
excitations is still of considerable current interest. 

Here we are interested in the Raman scattering from 
excitations of a 2-d normal electron system, namely from 
electron- hole pairs and collective plasmon excitations. 
We show that Raman scattering cross section in 2-d sys- 
tems differs from the spectrum of a 3-d metal in two 
aspects. First, the scattering from electron-hole pairs 
becomes more singular (due to the square-root singular- 
ity in the density of states). A finite strength of the 
electron-hole contribution is determined by the electron 
scattering rate. Second, as soon as the plasmon spectrum 
is gapless (uj oc k 1 / 2 ), a corresponding peak is located in 
the reasonably low-frequency (~ 10 meV) range. 

Inelastic light scattering in 2-d systems has been ex- 
tensively used in investigations of the excitations of the 
Fermi searJ )n the energy gappin the fractional quantum 
Hall regimea, excitpn stateso, spin-density and charge- 
density excitations! Experimental evidence of 2-d plas- 
mons in GaAs heterostructures comes from the Raman 
spectroscopy measurements in a magnetic fieldo. By 
varying the direction of magnetic field it is possible to 
distinguish Raman response of 3-2-d electron system from 
a contribution o£ a backgroundBl. The fc 1 / 2 -spectrum is 
clearly observedB, however, no data are available for the 
dependence of a lineshape of the plasmon peak on the 
electron scattering rate. 



The standard quantum mechanical theory of Raman 
scattering in electron systems applies the Green func- 
tion formalism. We use in the present paper a differ- 
ent approach, based on the Bqltpmann equation. Such 
a semiclassical method is validtil when the characteris- 
tic scale of transferred light momentum is less than the 
Fermi momentum of a 2-d electron gas. For the typical 
situation of GaA_/AlGaAs heterojunction the concentra- 
tion of carriersou is of the order of 1 — 6 x 10 11 cm -2 
and corresponds to the Fermi momentum value of ~ 
3 — 8 x 10 5 cm . On the other hand the typical val- 
ues of the momentum transfer are considerably lower 
~ 0.2 — 1 x 10 cm . The main advantage of the ki- 
netic approach is a possibility to include effects of elec- 
tron scattering and to account for electromagnetic field 
in a self-consistent manner by solving simultaneously the 
Maxwell equation. This gives a plasmon contribution to- 
gether with the electron-hole continuum. It is also possi- 
ble to generalise easily this theory to the case of external 
magnetic field applied. 
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FIG. 1. Geometry of the light scattering. Thick solid line 
(z — plane) represents the electron gas embedded into a host 
material with the dielectric constant e. The incident (i) light 
induces the in-plane 2-d current which in its turn produce a 
scattered (s) wave with different frequency. 

The system under study is shown in Fig. 1. Two- 
dimensional electron gas [z = plane) is embedded into 
a host material with the dielectric constant e. The po- 
larization of incident (i) and, hence, of scattered (s) light 
waves are assumed to be parallel to the plane z = 0. The 
effective Hamiltonian describing Raman scattering from 
electronic fluctuations is bilinear in the vector potential 
of light, 



H, 



eff 



2mc 2 



J d 2 s 5n 7 (s, t)A 2 (s, t), 



(1) 



where the fluctuation Sn-y expressed via the nonequilib- 
rium partition function Sf p (s,t) 



<5n 7 (s, i) 



2d 2 p 



(2) 



differs from the usual electronic density only by the 
anisotropic dimensionless factor 7 P (electron-light ver- 
tex). This factor depends on the light polarization and 
accounts for the virtual interband transitions. |J.ts exact 
form is not essential for the following (see Ref.E-3). 

Varying the expression ([l]) over the vector potential 
A we obtain the electron current induced by the incident 
light with the frequency a;W and the in-plane wave vector 

p 2 

j (i) (M) = dnJs,t)A^exp(-iw^t + ik^s), (3) 

mc 

The 2-d current (^) produces a scattered electromagnetic 

wave with different frequency a/ s ) and wave vector k s . 
The solution of corresponding nonuniform Maxwell equa- 
tion is straightforward. After some simple calculations it 
gives for the amplitude of light scattered into the half- 
space z > the expression, 



A(a,W J k«) = H^^^( £l ;,k.), 



mc 2 k 



(4) 



where ki = eu^'^/c 2 — kg ; the Fourier-component of 
the density fluctuations depends on the transferred en- 
ergy to — — lo^ and momentum k s = kg — kg . 

The Raman scattering cross section denned as the nor- 
malised energy flow ((\lu^ / lo^ | 2 )) related to 

the interval doj^ d 2 ki^ / (2ir) 3 , has the form: 



(Pa 



a/2,4 



,(s)3 



du^do^ 27rm 2 c 5 w (i)2^( s ) 



K(u, k s 



(5) 



where K(u>,k s ) is the Fourier component of the correla- 
tor of density fluctuations 

K(s -s',t- t') = {{5n*y(s, t)dn 7 (s', if))) . 



One can argue that the expression (0) diverges as the di- 
rection of scattered light approaches the electron plane: 

(s) 

fcz — > 0. In fact, this means that as soon as the "width" 
of a two-dimensional system I is assumed to be the small- 
est one of all of the characteristic lengths of the problem, 

(s) _ 1 

we are restricted to the limit k\ ' >> Z . 

To evaluate this correlator we apply the fluctuation- 
dissipation theorem which expresses it via the imaginary 
part of the generalized response Sn-y(uj, k) to an arbitrary 
external potential U(u>,k) (in what follows we omit the 
subscript s): 



K(tu,k) = - 



( <5n 7 (w, k) 
1 - exp (-w/T) V U(uj,k) 



Im 



(6) 



The most simple way to derive the generalised response 
(0) is to make use of the linearized Boltzmann equation 
for the nonequi] i br i urn part of the distribution function 
Sf P = x P dfo/d£M: 



— i(u/ — kv + ir 1 )x P (k,uj) = iujj p U(k,uj) 



evE(k,u;), 
(7) 



where fo(e) is the local-equilibrium Fermi-Dirac partition 
function. The second term in the right hand side of Eq. 
(0) accounts for the fluctuating electromagnetic field E. 
It satisfies the Maxwell equation with the nonequilibrium 
electric current determined from Eq. (0): 



rot rot E(z, s, w) — — 5-Efz, s, uo) 



(8) 



Here the brackets denote the integral over the Fermi line 



(...) = 



2dpp 
v(2ir) 2 



(...). 



We are interested in the solution of Eq. (g) at z = 0. 
The straightforward derivation gives for the Fourier com- 
ponent of the electric field 



E(z = 0,k,w) 



2ixie 



^/k 2 -eu 2 /c 2 <v Xp (k,w)>. (9) 



Substituting the solution (^) into the Boltzmann equa- 
tion ([7]) one gets the integral equation for the electronic 
density fluctuation x p (k, oj). Such an equation has a sim- 
ple solution which after the substitution into Eq. (|^) and 
then into the fluctuation-dissipation theorem @ gives 
the Raman cross section. Finally we obtain (see Fig. 2) 



K(k, oj) oc — Im 



U1 P 



, u> — kv + it 1 
Im F a (k,cj)D a0 {k,uj)F (k,uj) 



(10) 
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where the proportionality coefficient (Bose factor) is 
omitted, see Eq. (H); D(k, uj) is the two-dimensional elec- 
tromagnetic Green function 



2^e 2 



V a V/3 



uj \ uj — kv + IT 



2ne 2 v / k 2 -euj 2 /c 2 
(11) 



and F a (k,uj) is the oscillator strength 



— kv + IT 



(12) 






We devote the rest of the paper to the discussion of dif- 
ferent terms in the Raman cross section (llOf). 
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FIG. 2. Diagrammatic representation of the den- 
sity-density correlation function (Jioj) . The solid lines rep- 
resent electron propagators, the dashed line stands for the 
electromagnetic Green function and the empty vertex for the 
effective electron-light interaction. 

The first term represents the scattering from the 
electron-hole pairs. For the estimates we suppose the 
Fermi line to be isotropic, i.e., 



K e - h (k,uj) oc 



T«7 



Im 



y/(uj + IT- 1 ) 2 - k 2 V 2 



In the dirty limit kvr << 1 (the so-called 
momentum transfer limit) one can neglect v 



(13) 

zero- 
in the 



denominator. The cross section then.£akes the same well- 
known Lorentzian form as in 3-d caseO: ~ wt/(w 2 t 2 + 1). 

In the clean limit kvr >> 1 the expression ([l3]) has 
a square root singularity at uj — kv rather than a step- 
like one (as in 3-d case). It results in the strong non- 
symmetric resonance (see Fig. 3); the finite height of 
this resonance is controlled by the scattering rate t . 
For anisotropic Fermi line the resonance location is de- 
fined by the maximum value of electron velocity along 
the momentum transfer u> = ky mn3: . 

The second term in Eq. ([H]) represents effects of 
Coulomb interaction and collective electron excitations, 
namely 2-d plasmons. This term is important in the clean 
limit only. At low transferred frequencies uj << kv it re- 
sults in the screening of the isotropic scattering channel 
through the renormalization of the electron- light- vertex 
7p - ► 7p — (7p) / (1)) similar to usual 3-d caseOEj. 

At high transferred frequencies uj >> kv the second 
term in Eq. (|o|) gives a 2-d plasmon peak located at the 
plasmon frequency u> p i(k) which is determined by the dis- 
persion equation 



(v 2 k ) sjk 2 - euj 2 /c 2 , 



(14) 



where Vk means the component of electron velocity in 
the /c-direction. The typical momenta transfer k are of 
the order of light momenta. Hence from the formula Jl4| ) 
one can see that u> » vk proving the initial assump- 
tion was correct. We can also omit the term euj 2 /c 2 
(this term cares for the finite plasmon velocity) in Eq. 
(n~4|) in comparison with the term k 2 due to the fact that 
<rk >> v 2 pf for typical values of k. Indeed, this means 
that it is enough to use the Poisson equation for electro- 
magnetic fluctuations instead of the Maxwell equation 
(g). The only difference occurs at very small transferred 
momenta where the Poisson equation gives the infinite 
plasmon velocity in the limit k — ► 0. 

The formula (|l4|) is valid only if the plasmon wave- 
length becomes large compared to the layer thickness, 
k >> If this condition is violated the three- 

dimensional problem has to be solved with boundary con- 
ditions satisfied on both sides of the layer. It's solution 
gives the expression 



4») 



where the electron velocity along the perpendicular di- 
rection v z appeares. Note, that the angular brackets now 
denote the integral over the three-dimensional Fermi sur- 
face. When I — y this expression reduces to Eq. ( |l4| ) and 
for / — > oo it gives the frequency of ordinary 3-d plasmon. 
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FIG. 3. Raman scattering cross section. The 

non-symmetric resonance is located at uj — kv. Its width at 
uj > kv as well as its height are determined by the finite value 
of the scattering rate t~ . The second symmetric resonance 
corresponds to excitation of plasmon. 

Near the plasmon resonance the Raman cross section 
has the symmetric Lorentzian lineshape 



K p i(k,uj) oc 



7717 



StTCJ (Wpi(fc) - uj) 2 +t- 2 /4' 



(15) 
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The relative height of two resonances ( |l3| ) and ( |i~5| ) is 
Ke-h/Kpi ~ fc 3 / 2 t; 3 / 2 r 1 / 2 /a;pj and can be either more or 
less than unity depending on the momentum transfer k 
and scattering rate r _1 . 

In conclusion, we have calculated the Raman scat- 
tering intensity from two-dimensional electronic fluctu- 
ations. The main distinctive features from a usual three- 
dimensional metal are: the more singular electron-hole 
contribution ( |l3| ) and low frequency plasmon resonance 
( |l5| ) . The electronic Raman scattering in a 2-d system in 
a transverse magnetic field can be considered studied the 
same Boltzmann equation technique as it has been done 
for 3-d electron systemEl 
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